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Abstract 

In relativistic heavy ion collisions, also the bremsstrahlung of gluons in the frag- 
mentation regions of the nuclei suppresses the produced charmonium states. In the 
energy range of the SPS, the radiation of semi-hard gluons occurs in the Bethe-Heitler 
regime and the density of gluons and therefore the suppression goes like (AB) 3 , where 
A and B are the nucleon numbers of the projectile and target nuclei. In contrast, 
the suppression via collisions with nucleons is proportional to (As Parameter 
free perturbative QCD calculations are in a good agreement with the data on J/^f 
and suppression in heavy ion collisions at SPS CERN. At higher energies (RHIC, 
LHC) the number of gluons which are able to break-up the charmonium substantially 
decreases and the additional suppression is expected to vanish. 



Challenges in charmonium production off nuclei 

Charmonium production in heavy ion collisions is considered as one of several sensitive 
probes for the formation of a quark-gluon plasma |l] . This idea has led the NA38 and NA50 
collaborations at SPS CERN to perform a series of measurements of nuclear suppression for 
charmonium production in p A and A B collisions. We recall some remarkable observations 
(see the collection of data in 0): 

(i) If J/ty suppression in proton-nucleus collisions is treated as a final state absorption 
phenomenon one extracts an absorption cross section <j in (J /'ty N) k, 6 — 7mb from the data. 
At the same time, the vector dominance model applied to J/ty photoproduction data on 
protons results in a much smaller value cr to t(J r /\I/ N) ~ 1.1 mb, which increases to a value 
of 3 — Amb for a coupled-channel analysis of the photoproduction data ||. However, a 
discrepancy of about 3 mb still remains to be explained. A part of this extra suppression of 
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J/^f may arise due to a stronger absorption of Xc states which contribute to the production 
of J/ty about 30% via decayQ. 

(ii) The ^'(2S) state is expected to have a radius about twice as large as the J/\l/ and, 
therefore, to have a few times stronger final state interaction. However, in proton nucleus 
collisions both J/\l/ and ^' are suppressed by nearly the same amount. This is an effect 
of the finite formation time of the J/\l/ and \P' and can be quantitatively explained in the 
coupled channel approach ||. 

(iii) However, in S-U and Pb-Pb collisions ^' is much stronger suppressed than J/^>. 
This is partially an effect of the formation time together with the inverse kinematics for the 
projectile nucleus |§. Nevertheless nearly a half of the observed suppression still needs to 
be explained. 

(iv) Final state absorption on nucleons well describes the data on J/^ suppression in 
p-A, O-Cu, O-U and S-U collisions with light projectiles. However, it fails in the case of 
Pb-Pb collision where data show a much stronger suppression of J/^f. This observation 
presented at QM'96 [[| and published in has led to intensive theoretical work (for a 
recent review we refer to [0]). At present there are two main schools of thought: A number 
of authors consider the strong suppression in the Pb-Pb data as a proof for the formation of 
a QGP, in the spirit of the original idea by Matsui and Satz M. Other authors explain the 
strong suppression of J/ty and *f>' in Pb-Pb by hadronic "comovers", i.e. via destruction of 
the charmonium by the hadrons which are produced in the final state of a nucleus-nucleus 
collision with the same velocity as the charmonium. Both approaches, via QGP or comovers, 
introduce a number of free adjustable parameters, and it is only then that they can account 
for the data. For this reason it is not possible yet to decide about the mechanism. 

In this paper we are able to explain the same set of data by invoking - an addition to the 
normal suppression via collisions with nucleons - the charmonium suppression due to gluon 
bremsstrahlung produced in multiple nucleon interaction. Since we deal with semi-hard 
processes, pQCD is expected to work reasonably well and no adjustable parameters have to 
be introduced. 

Charmonium break-up via interaction with nucleons and gluons 

Final state interaction of a charmonium in nuclear medium which leads to the nuclear 
suppression is usually described in terms of collisions of the charmonium with undisturbed 
bound nucleons. This is might not be fully correct. A nucleus-nucleus collision is illustrated 
in Fig. |l| on a two dimensional (time - longitudinal coordinate) plot which corresponds to the 
NN cm. frame. One can see that some of the nucleon trajectories cross the charmonium 
after they have interacted with nucleons from another nucleus. A natural question arises, 
whether the charmonium interacts with such debris of a nucleon in the same way as with 
an intact one? 

To answer this question note, that a small-size cc pair interacts mostly with only one 
of the valence quarks in a large-size proton. This is because of color screening which cuts 
off the gluons propagating far away from the charmonium (the Van der Waals forces are 
supposed to be cut off by confinement). Thus, the charmonium acts like a counter of the 
number of constituents in the nucleon. In the constituent quark model we expect 

<r in (*N) «3<7 in (ttg) . (1) 
1 We are thankful to Mark Strikman for this comment 



2 




B A 

> z 

Figure 1: A two-dimensional (time - longitudinal coordinate) plot for 
charmonium production in a collision of nuclei A and B in the cm. of 
the colliding nucleons. The solid and dashed lines show the nucleon and 
gluon trajectories, respectively. 



Interacting with the debris of a "wounded" nucleon (N*) the charmonium can find more 
constituents in it than in the original nucleon, particularly gluons shown by dashed lines 
in Fig. [I], which had a sufficiently short radiation time. In this case the cc break-up cross 
section increases, 

a in (*N*) = a in (^N) + a in (%) (n g ) , (2) 

where (n g ) is the mean number of gluons radiated in a NN interaction preceding the collision 
with the charmonium. 

Interaction of the cc pair with "wounded" and intact nucleons leads to the following 
suppression of charmonium in AB collision, 
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Here \I/ denotes either J/ty or Effective absorption cross sections of the cc pair a e jf(^N) 
and a e ff(^fg) may be different for nuclei A and B (for nonzero and depend on where 
and when the J/\I / or ty' is produced (see below). 

In Eq. (H) the coordinates of the production point for the cc pair are assumed to be 
(s, za) and {b — s, zb) in the nuclei A and B, respectively, and b is the impact parameter for 
the AB collision; pa,b{^} is the density of the nuclei A or B. 

The first exponential in ([|) accounts for a suppression due to interaction of the char- 
monium with the nucleons, both in nuclei A and B. The second and the third exponential 
factors include suppression of the charmonium due to interaction with the gluons radiated 
in the fragmentation regions of the nuclei A and B, respectively. The last exponential elim- 
inates a double counting in number of interacting nucleons (assuming for this purpose that 
the charmonium is produced with xp — 0). 

Gluon bremsstrahlung will also affect nuclear suppression in proton-nucleus collision. 
Reducing Fig. 1 to the case of pA collision one can see that some of gluons radiated by 
the projectile proton cross the trajectory of the charmonium. The expression for nuclear 
suppression is simpler than @, 

^ oo oo 

S^ A = — j d 2 b J dz A Pa(P, z a ) exp -af ff (^N) J dz p A (b, z 
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Here we neglect the possible effect of dilution of the gluon density due to transverse motion 
of the gluons, since for the time interval less than one Fermi most of them stay within the 
nucleon size in transverse plane. 

In order to evaluate expressions (§) - ©, we have to discuss the input parameters, the 
mean number of produced gluons (n g ), the charmonium break-up cross sections <Ji n (^N) 
and <Ji n (tyg) which are effective cross sections because of the formation time effects. All the 
quantities can be estimated fairly reliably. We concentrate on the kinematics of SPS which 
corresponds to about 10 GeV per nucleon in cm. frame. All the longitudinal distances in 
the nuclei are contracted by a factor of 7 = 10. 

Mean number (n g ) of radiated gluons 

A quark in a projectile nucleon radiates only if interacts with a target nucleon ||. The 
mean free path for a quark in nuclear medium is three times longer than for a nucleon, 
A 9 ~ 6 fm. One has a maximal number of produced gluons, namely, the transversal density 
of gluons proportional to (AB)s, if the gluons are radiated independently in each NN 
interaction in Fig. [l|. This is possible only if the radiation happens in the Bethe-Heitler 
regime, i.e. if the formation length l 9 j of the radiation in the cm. frame does not exceed 
the mean free path of a quark Az = Ag/7 ~ 0.6 fm: 

^ a i1':°' ^, <5) 



a 2 mj, + k 2 



<2 



where E q and m q are the initial energy and the mass of the radiating quark; a and k are 
the fraction of the quark light-cone momentum and the transverse momentum carried by 
the gluon. 
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When the condition (|5]) is violated, the Landau-Pomeranchuk effect suppresses gluon 
radiation, besides, most of gluons are formed too late to cross the charmonium trajectory 
in Fig. [l|. 

Thus, the mean number of gluons radiated in a NN collision is given by the integral 
over the gluon radiation spectrum weighted by a step-function which requires the radiation 
to be in the Bethe-Heitler regime, 
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Here a min = (u min + \ju 2 min - k 2 )/(2E q ) if k < u min , otherwise a min = k/(2E q ). The 
choice of the soft limit k min has practically no influence on (n g ) because the step function 
in (|6|) cuts off the gluons with small transverse momenta. Those which contribute have 
k > 0.3 — 0.5 GeV. These are semi-hard gluons which do not overlap with the soft gluonic 
field in a color-flux tube (string). 

The cross section of gluon radiation as function of a and k integrated over the final 
transverse momentum of the quark is calculated in U, 



da(qN gX) 3 a s (k 2 ) C 2 m\ a 4 k 2 + [1 + (1 - a) 2 } (k 4 + a 4 m 4 ) 
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Here a s (k 2 ) is the QCD running coupling; C is the factor for the dipole approximation for 
the cross section of a qq pair with a nucleon, <7 59 (r<r) ~ C rf., where Tt is the qq transverse 
separation |I| ||. PQCD predicts C ~ 3. Note that (0) is supposed to work well even for 



radiation of gluons with rather small k 0. 

We are left with only one parameter io min in @ which brings the main uncertainty in 
the value of (n g ). In accordance with the consideration in next section we try two values 
Umin = 0.5 and 1 GeV, which result in (n g ) = 0.69 and 0.25, respectively. 

Note that the radiation spectrum steeply grows with energy since we select gluons with 
relatively large kr- According to the analysis [[11], |9| of HERA data for the proton structure 



function we expect energy dependence dn g /dadk oc (s/s ) . However, the restriction for 
the formation length imposed by the step function in (|6|) suppresses (n g ) much more at high 
energy because of Lorentz contraction of the nuclei. 

We compare (n g ) for AB collision calculated for uj m i n = 0.5 GeV at SPS, RHIC and 
LHC energies including also the energy growth of the gluon spectrum. 

f 6.9 x 10- 1 (SPS, v/i = 20 GeV) 
(n g )= I 6.9 xlO- 3 (RHIC, y/s = 200 GeV) . (8) 
[l.2xl0~ 3 (LHC, y/s = 1200 GeV) 

In the rest frame of either of the colliding nuclei the main fraction of gluons are radiated at 
high energies later than the charmonium is produced. 

The inelastic gluon - charmonium cross section (Ti n (^g) 

The choice of u min deserves a special discussion since it is related to the problem of 
the energy dependence of J/ty — g break-up cross section a. in (J /^). At high energies it is 
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dominated by gluonic exchange in t-channel and slightly grows with energy ||. According 
to the previous discussion and Eq. (|1|) the inelastic cross section with a quark is known. 
The one with a gluon differs by a color factor 9/4. Therefore, we expect 

a in (J/^g)^^a in (J/^N) . (9) 

In the energy range of interest a coupled-channel analysis of J/ty photoproduction data 
gives cr in (J/^/N) « Amb. Correspondingly, a in (J/^/g) ~ 3mb. For ip'{2S) we expect a 
larger cross section. If the cross section is proportional to (r 2 ) of the charmonium state, one 
has a factor of 7/3 in an oscillator model, or a factor of 4 for (r 2 ) calculated from realistic 
wave functions |12| . 



One could expect a vanishing &i n {tyN) down to the threshold energy which is about 
0.6 GeV lower for than for J/^J. However, since radiation time is of the order of nuclear 
radius, the gluons cross the charmonium trajectory in Fig. 1 at about the same time as the 
wounded nucleons, i.e., when the charmonium is not yet formed. One cannot require energy 
conservation to hold with accuracy better than AE w 1/Af where At is the time interval 
between the production moment of cc (the production/coherence time is quite short at the 
SPS energy) and its crossing with a gluon. The mean time of crossing in Fig. 1 is At « 
(1/3)Ra/j ~ lGeV" 1 . Therefore, one can disregard energy conservation, binding energy 
and the difference in threshold energy between J/ty and ty' with accuracy AE ~ lGeV. 
For the same reason, the energy threshold has no influence on ai n (^ g). 

There is also another mechanism of charmonium break up by gluons via direct absorption 



of the gluon which has a pick about ai n (^> g) ~ 5mb near the threshold with a width 
of the order of binding energy and falls steeply at higher energies. In this case one should 
replace the binding energy by the AE. 

To simplify calculations we fix <r in (J/^g) rs 3mb at uj > 0J min . We try uj min = 0.5 and 
1 GeV in accordance with the scale imposed by the estimated uncertainty in energy. 

Effective absorption cross sections 

The produced cc colorless pair propagates through the nucleus with a smaller size and 
a smaller absorption than that for either of J/ty and unless time exceeds the formation 
length Ij for the charmonium wave function |TJ], || , 

7* - 2E * nry) 
l f ~ Ml - Mj/„ • (10) 

Here Ey is the charmonium energy in the rest frame of the nucleus. The transition between 
short and long times is controlled by the nuclear formfactor |4j], which for the charmonium 
production point with impact parameter s reads, 



F A {q L ,s) = —- dze^ z p A (s,z) , (11) 

where the longitudinal momentum transfer is related to dlPf) , qi = 1/^?, and T(b) = 
J^oo dz p A (b, z) is the nuclear thickness function. 
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The break-up of the colorless cc pair during its evolution is described by the effective 
cross section M, 



a 



tf f {J/*N) = a in (J/^N)[l + eRF'i(qt,s)\ (12) 



atjjV&'N) = a in (^N)[l + —Fl(qi,s)\. (13) 
The parameters here are defined in ||. We use the values given by the harmonic oscillator 



model, r = (y'\&\y')/(J/y\a\J/V) = 7/3, e = (¥\a\J/^f)/(J/^\a\J/^) = -^2/3. The 
relative production rate of ty' to J/ty in a NN collision is known from experiment, R ~ 0.4. 

We expect different effective cross sections for absorption in the nuclei A and B, since 
the charmonium with (xp) = 0.15 (NA38, NA50) propagates through these nuclei with 
different energies in their rest frames, i.e. and q^ are different. In an AB collision with 
Ei a b = 158 GeV / A&nd Xp = 0.15 one has if — 2 fm and lp — 4 fm in the rest system of the 
projectile (A), target (B), respectively. The values for the effective cross sections in a Pb-Pb 
collision are (Eqs. (03)- QI^)) a£ ff (J/^N) = 3.9 mb, af ff (J/^N) = 3.3 mb, <rf ff (V'N) = 
8.6 mb and a? ff (WN) = 4.8 mb (for a in (J/^N) = Arab and a in (^'N) = 9.3mb) 

Due to the uncertainty for the radiation time one can assume that the gluons radiated 
by the "wounded" nucleons cross the $ trajectory at about the same points as the nucleons. 
Therefore, we can use the same form of the effective cross sections flT2] ) - ([ll$p for interaction 
with gluons. 

Results of calculations 

To simplify the calculations we use a constant nuclear density p(r) = p ®(Ra — r ) with 
Po = 0.16 fm~ 3 , except the nuclear formfactor where a realistic density is used. We calculate 
(0) and (U) using the values of the parameters fixed in previous sections. We use two values 
of (rig) = 0.25 and 0.5 to characterize of the uncertainty in our calculations. 

The results are plotted in Fig. 2a for J/ty and Fig. 2b for as function of A x B 
together with available experimental data for pA and AB collisions. We reproduce the 
data for J/ty suppression quite well, including the S-U and Pb-Pb points. We obtained a 
reasonable agreement with the data for suppression as well, although the curves seem to 
be a bit too high. In this respect we should note that the two coupled channel model we 
use with oscillator potential might be a poor approximation in the case of 1 i r '. Indeed, as 
was mentioned, the ratio of mean square radiuses of ^' to J/^f which is 7/3 in the oscillator 
model, is predicted to be nearly 4 in a more realistic approach. As an example we tried 
another value r = 3 in (|i~2f) - (|13D . The result is shown by the dashed curve in Fig. 2b which 
is in a better agreement with the data. Our results for J/ty suppression are very stable 
against this modification. 

Note that our curves are not straight lines even for pA collisions. This is mostly due to 
the formation time effects, i.e. is a result of A-dependence of the effective cross sections 
dT2] ) - (|T^). Furthermore, we are able to account for the J/*f> suppression in pA collisions, 
although we use <Ji n (J/tyN) = 4mb. The additional suppression (even in pA) comes from 
the gluons with absorption cross section of 3 mb. At larger values Ax B for nuclear collisions 
our curves deviate even more from straight lines (in agreement with the data!) due to the 
nonlinear dependence on nuclear radius, as is brought into by the gluon radiation, whose 
density is proportional to (AB) 1 / 3 . 
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Figure 2: (a) Nuclear suppression of J/ty production in pA and AB 
collisions as function of the product A x B. The two curves correspond 
to (n g ) = 0.25 and 0.5 The lines connect the point calculated for A x 
B corresponding to each experimental point, (b) The same as on (a), 
but for production rescaled by factor 10~ 2 . The circles and squares 
correspond to pA and AB data, respectively. The data points are from 

g H W- 



Formulas (^) and can be also used for calculation of nuclear suppression of charmonia 
as a function of impact parameter or of the mean length L of the total path of \1/ in nuclear 
matter. Our results for S-U and Pb-Pb collisions are plottes in Fig. |3] by dotted and solid 
curves respectively. We use (n g ) = 0.5 at 200 GeV and (n g ) = 0.57 at 158 GeV in accordance 
with Eq. (|^). The overall normalization is free. 

Following recommendation of the referee for our paper we compare the calculations with 
the published data [|J depicted in Fig. |] by round and triangle points for S-U and Pb- 
Pb collisions respectively. However, the assignment for L corresponding to the measured 
transverse energy is not trustable. It is based on an oversimplified model, particularly 
neglecting the Landau-Pomeranchuk suppression (at low k) and enhancement (at high k) 
of the gluon radiation spectrum |J. The cascading of gluons radiated in the Bethe-Heitler 
regime which increases the transverse energy, is also neglected. Besides, the normalization 
on the cross section of the Drell-Yan reaction brings an additional uncertainty (especially 
at high energies). We plan to do our own calculations for the transverse energy produced 
in heavy ion collisions which would allow an unbiased comparison with data. 

Conclusions 

Charmonium suppression in pA and AB collisions at SPS energies can be calculated as 
arising from two sources: 
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Figure 3: Nuclear suppression of J/^ relative to Drell-Yan lepton pairs 
as function of the mean length of total path of the charmonium in the 
colliding nuclei. The dashed (round points) and solid (triangles) curves 
correspond to S-U and Pb-Pb collisions, respectively. The data points 
are from ffij. 



(i) Collisions with nucleons, the cross section for which is the one deduced from photo- 
production data modified because of the formation time effects. The suppression from this 
source essentially scales like A 1 / 3 + B 1 / 3 . 

(ii) Collisions of charmonia with gluons which have been produced in multiple NN 
collisions. In the Bethe-Heitler regime - good at SPS energies - about (n g ) = 0.25 — 0.5 
gluons contribute per NN collision. This additional effect scales with (AB) 1 ^ 3 and accounts 
for the "anomalous" J/\I/ and suppression in nucleus-nucleus collisions. We conclude 
that the results of the NA38/50 experiments do not signal about QGP formation. 

Nonetheless, the suggested mechanism leaves room for other contributions to the char- 
monium suppression within the uncertainty of calculations. Particularly, as was mentioned 
above, a "wounded" nucleon is not a colorless system of partons any more. Therefore, color 
screening does not cut off (only the charmonium formfactor does) soft gluon exchanges in 
contrast to interaction. However this correction to ([J), which we do not expect to be 
large, scales with {A 1 / 3 + B l l 3 ) and leads to a renormalization of ai n {9N). 

According to (||) the suppression of charmonium by gluon radiation vanishes at high 
energies due to the formation time effect. We expect practically no additional suppression 
by prompt gluons at RHIC, while the suppression via interaction with "comovers" (either 
QGP or hadrons) increases with energy. Therefore, one has less freedom in interpretation 
of data on nuclear suppression of charmonium at RHIC. 

The main parameter which controls the charmonium suppression via interaction with 
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prompt gluons in (|) is the product (n g )a e ff(^N). We estimated this factor and allowed it 
to vary within a factor of 2. We think that this uncertainty covers other possible corrections, 
for instance, the contribution of Xc decays to the production of J/^f. 
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